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Abstract 

We discuss the constraints induced by the algebra of the Poincare generators on 
non-relativistic effective field theories. In the first part we derive some relations 
among the matching coefficients of the HQET (and NRQCD), which have been 
formerly obtained by use of reparametrization invariance. In the second part we 
obtain new constraints on the matching coefficients of pNRQCD. 
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1 Introduction 



For any Poincare invariant theory the generators H, P, J, K of time translations, space 
translations, rotations, and Lorentz transformations satisfy the Poincare algebra: 



[P\ p j ] 


= o, 


(1) 


[P\H] 


= o, 


(2) 


[J\ P J ] 




(3) 


[3\H] 


= o, 


(4) 






(5) 


[P\ K j ] 


= -iSijH, 


(6) 


[H,K l ] 


= -iP\ 


(7) 


[J\K j } 




(8) 


[K\ K j ] 


-ik 


(9) 



It has been pointed out, as early as in Ref. [1], that the algebra induces non trivial 
constraints on the form of the Hamiltonian of non-relativistic systems where Poincare 
invariance is no longer explicit. Indeed, the algebra has been used in the past to constrain 
the form of the relativistic corrections to phenomenological potentials [2]. 

In this letter by Poincare invariance we mean the explicit realization of the algebra 
(l)-(9). We study such realization in some of the modern non-relativistic effective field 
theories of QCD: Heavy Quark Effective Theory (HQET) [3], Non Relativistic QCD 
(NRQCD) [4] and potential NRQCD (pNRQCD) [5, 6]. One may expect that also 
in these cases Poincare invariance induces non trivial constraints on the form of the 
interaction. More specifically, one expects to obtain some exact relations among the 
matching coefficients of the considered effective field theory. This study has never been 
done before. In the case of the HQET some exact relations between matching coefficients 
have been derived using a specific invariance of the theory, known as reparametrization 
invariance [7, 8]. We will derive in Sec. 2 some of these relations, showing in this way 
that reparametrization invariance is just a manifestation of the Poincare invariance of 
the theory. In the case where dynamical gauge fields have been integrated out from 
pNRQCD, Poincare invariance has been studied in this framework in [9]. There some 
relations among the potentials have been obtained that were previously known only from 
the transformation properties under Lorentz boost of the representation of the potentials 
in terms of Wilson loops [10, 11, 12]. Here we will extend that study to the situation 
with dynamical gluons. This is the situation of interest, for instance for weakly coupled 
heavy quarkonium states like the T(IS'). 

The letter is organized as follows. In Sec. 2 we illustrate how Poincare invariance 
(in the sense specified above) works in NRQCD/HQET. In particular we construct all 
generators of the Poincare transformation up to the relevant order and verify the algebra. 
We obtain in this way some of the constraints already known from reparametrization 
invariance. In Sec. 3 we apply the same machinery to pNRQCD. After constructing 
all generators of the Poincare transformations up to the relevant order, we derive, by 
imposing the algebra, some new constraints on the matching coefficients of the theory. 
In Sec. 4 we summarize our results and discuss some possible future developments. 



1 



2 NRQCD 



After integrating out the hard scale m from QCD, one obtains NRQCD [4]. Neglecting 
operators that involve light-quark fields [13], the most general NRQCD Lagrangian den- 
sity (up to field redefinitions) for a quark and an antiquark of mass m up to order 1/m 2 
is given by (we display also the term D 4 /(8m 3 ) for further use): 

£nrqcd = 

D 2 D 4 (t-B . [D-,E] <t-[Dx,E]1, 

^Y D °- m + Cl ^ + C2 8^ + CF9 ^ + CD9 ^^ + tCs9 8m 2 f 
, i D 2 D 4 cr B . [D-,E] «t-[Dx,E] 

+X \ tDo + m - Cl 2^- C *8^- CF9 ^ + Cd9 ^T + lCs9 8m 2 J* 

+—^Wx + ■ X^X + ^ f T>x f T a x + %V f T a ^ • X f TV X 

m z m z m z m z 

--F; v F a ^ + %f abc F; v Fl a F^ (10) 

t: lit 

where ip is the Pauli spinor field that annihilates the fermion and \ is the Pauli spinor field 
that creates the antifermion, iD° = id°-gA°, iD = iV+gA, [D-, E] = D E-E • D, and 
[Dx, E] = D x E - E x D, E* = F i0 and B i = -e ijk F jk /2, e ijk being the usual three- 
dimensional antisymmetric tensor ((ax b)* = e^a 7 !^). We display the mass terms, that 
are usually removed by a field redefinition, for the further use of the canonical formalism 
to derive the time-translation generator. The one-loop expressions in the MS scheme for 
the coefficients cp, c' D , c$ and d' 3 can be found in [8] (according to the definitions of c' D 
and d' 3 given in [14]) and for dij = s,v) in [15]. 

It has been proved in [8] that in the bilinear sector the NRQCD Lagrangian is equiv- 
alent to the Lagrangian of the HQET. Therefore, at the order at which only the bilinear 
sector matters (and this will be the case for the rest of this section), the results that 
we obtain are valid both for NRQCD and HQET. The main reason for introducing here 
the NRQCD Lagrangian up to order 1/m 2 , is that it is needed to obtain the pNRQCD 
Lagrangian at 0(l/m 2 ), which is the subject of Sec. 3. 



2.1 Canonical quantization 

To be definite we quantize NRQCD in the A = gauge [16]. The pairs of canonical 
variables are (x^X^) an d (Aj a , II^ = <9£nrqcd/ d(d A^)). The physical states 



|phys) are constrained by the Gauss law: 

(D • nr|phys> = <?(^T> + X t T a x )|phys). (11) 
The canonical variables satisfy the usual equal time commutation relations: 

[Ili(x,t),Ai(y,t)] = ^.(U^(x-y), (12) 

[A£(x,t),Aj(y,f)] = [n*(x,t),nKy,t)] = o, (13) 

{^(x,t),^( y ,t)} = { XQ (x,t),xJ j (y,t)} = ^ (3) (x-y), (14) 

0Mx,t),^(y,t)} = {^(x,t),^(y,t)} = 0, (15) 

{Xa(x,t),xp(y,t)} = {xt(x,*),xj,(y,t)} = 0. (16) 
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All other commutators are zero. In the following, in order to fulfil the Gauss-law con- 
straint (11), we will assume the commutators to act on a space spanned by the physical 
states. 



2.2 Poincare algebra generators in NRQCD 

The construction of the generators proceeds in the following way. The generators H, 
P and J can be derived from the symmetric energy-momentum tensor (see for instance 
[17]). Since translational and rotational invariance remain exact symmetries when going 
to the effective theory, the transformation properties of the new particle fields under 
these symmetries are the same as in the original theory. The derivation of the Lorentz- 
boost generators is more problematic, since the non-relativistic expansion has destroyed 
the manifest covariance under boosts. A consistent way to construct K is to write down 
the most general expression, with some obvious restrictions already included, containing 
all operators consistent with its symmetries and to match it to the QCD Lorentz-boost 

r \ ( Yl a ■ IT a -|- B a • B a 1 

generator, which is K = — t P + J d 3 x - < x, h *0(— iD • 7 + m)ip >. 

This is very much the same procedure that is used in the construction of the NRQCD 
Lagrangian. Accordingly, new matching coefficients, typical of K, will appear. 

Specifically for the Lagrangian (10) we obtain (P and J are exact, H is displayed 
up to order 1/m 2 , including the kinetic energy up to order 1/m 3 , and K is given up to 
order 1/m): 

h = Am - Cl — - c 2 — - c F g— - c' g [B ' ,Il] - iz ° ' [D>< ' U] )j> 
\ 2m 8m 3 2m D 8m 2 8m 2 J 

J_ m + c ^! + c jE!1 + c a B c > [D-.n] ic o-- [Dx,n] \ 

X \ 1 2m 2 8m 3 F ^ 2m D ^ 8m 2 8m 2 J ^ 

-V^X - — 2 ^^ ■ X^X ~ —^ a ^ a X ~ — 2 ^*i> ■ X t T V X 
m A m z m A m z 

™ ™ , do . pa J/ 

I m E = II 

H = J d 3 xh, (17) 

P = Jd 3 x ^(_,D)^ + x t(_,D) x + i[n a x,B a ]), (18) 
J = J d 3 x (V (x x (-iD) +^ + X f (xx (-iD) + |) X 

+ixx[ITx,B a ]), (19) 
K = -tP + Jd 3 x^± 

-*« / d 3 x (±r> I x (-iD) , x (-iD) X ) , (20) 

where A;* 1 ) is a matching coefficient specific of K. 
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2.3 Poincare algebra constraints in NRQCD 



Let us now consider the constraints induced by the Poincare algebra (l)-(9) on the 
NRQCD generators H and K. The constraint [P\ K J 'l = —i5nH has been already 

used in Ec, (20). Indeed, this commutation Nation tees K to have the form j A 

{x, h(x, t)}/2 + translational-invariant terms that depend on x only through the canon- 
ical variables. From [K*,K J ] = —ieijkJ k at O{l/m ) it follows that 

-ie ijk (1 - feW) J d 3 x (V ^ + x^x^j = =► fc« = i. (21) 

From [H, K<] = -iP i at O(l/m ) it follows that 

-i (1 - Cl ) J d 3 x (f HD l ) i, + x t (_,d j ) X ) = Cl = 1, (22) 
and at 0(l/m) 

-i(2c F -c s -l) /A ^.^^V'-X t ^|^x)=0^2e F -e s -l = 0. (23) 

Finally from [H, K*] = -iP* at £>(V 2 V7m 2 ) we obtain 

(l-c 2 ) y^^t__^ + x t__ x j =0 =► c 2 = l. (24) 

All other commutation relations are satisfied at the order we are working. The constraints 
(22), (23) and (24) were first derived in the framework of reparametrization invariance 
in [7, 8]. 



3 pNRQCD 

The pNRQCD Lagrangian for a heavy quark-antiquark system is obtained from NRQCD 
by integrating out the soft degrees of freedom associated with the scale of the relative 
momentum of the two heavy quarks in the bound state [5, 6]. The name pNRQCD has 
been used in the literature to identify effective field theories with different degrees of free- 
dom. Here we call pNRQCD the effective field theory that can be obtained from NRQCD 
by perturbative matching and contains, as degrees of freedom, the quark-antiquark field 
(that can be split into a colour singlet S = S\ c /\fW c and a colour octet O = O a T a / ' \fTp 
component) and (ultrasoft) gluons. The fields S and O a are functions of (X, t) and x, 
where X = (xi + x 2 )/2 is the centre-of-mass coordinate and x = xi — x 2 the relative 
coordinate, with x x (x 2 ) the coordinates of the quark (antiquark). The coordinate x 
plays the role of a continuous parameter, which specifies different fields. All the gauge 
fields have been multipole expanded around the centre-of-mass and depend on (X, t): 
p^ua = i?^«(x,i) and iD^O = id^O + ig[A^(X, t), O]. The terms in the pNRQCD 
Lagrangian are organized by powers in the 1/m and x expansions. We will indicate a 
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generic term of order x 3 /m ! as where ip, <fi G {S,0}. Notice that the derivative V x 
counts like x~ x and, therefore, contributes with a negative power to the second index of 

n <p4> ■ 

We aim at verifying the Poincare algebra at order x/m° and x°/m. For this purpose 
we need the pNRQCD Lagrangian at order x 2 /m° (obtained in [5, 18]), as well as at order 
x°/m, (x/m) Px and (x°/m 2 ) Px (derived in this work). We include all local operators 
of the correct transformation properties and the appropriate dimensions, moreover we 
assume that, apart from the kinetic energy, momentum and spin operators are suppressed 
by a factor 1/m each. To the above order the pNRQCD Lagrangian density is given by 
(as in NRQCD, we display the mass terms, that could be removed by a field redefinition, 
for the further use of the canonical formalism to derive the time-translation generator): 

£ p nrqcd = / d 3 x Tr |s f (ido -2m-h s )S + O f (iD - 2m - h Q ) O 

- [(S f /i5oO + H.C.) + C.C.] - [0 ] h 00 + C.C.] - [0 ] h% 0h% + C.C.] j 

-\f; v F^\ (25) 



= Vrip {x) + 1 -{P\, V p ^ a (x)} + l -{ V l V p ^ b (x)} 

(x x P x ) 2 . (x x p x ) 2 
H — Vv^x) H = Vjjt^x) 



with 



(26) 



Ax 2 



x z 



+ (xxPx)-(^>-^) (g) + (xxp.)-(^) + ^) 
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x ■ er^ x • a 



ar 



h 



7,(0,1) ,(0,2) , .(1,0) ,(l,l)/ p x , 1,(2,0) / p \ 

x-#E, 



2 

(2) 

with 



<7« ■ <T< 2 > 



1.(0,1) 



-(0,1), 



(0,2) = 

<p<t> 
1/1,0) _ 



K0,2), 



- Va o l j xV(D'gE^) - ^ & 1 V (D • <?E), 
8 8 

^ Vfif (x) { P „ x x gB} - V^f>(x) • *B 



2m x 2 
1 



(x.^)(x.^B) 



m 2a; 



x-#E, 



Px) = t— KJ 1} (:t) {P r ,xx 5 B}, 



h 



(2,0), 



Py) = 



8m 

c s 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 



vS?(x)<rU.[P x x,gE} 



16m 2 
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1 V™{x) 



16m 2 x 2 
1 V*${x) 



(x^ (1) ){Pr,(jExx)} 
{(x^E),P x .(xxo-«)} 



16m 2 x 2 

^ W {(x-Px),^.(xx // E)} 



1 



16m 2 x 2 
T ^{(p,-Px),^ ) ( a ;)(x.,E)} 



16m 2 
1 



16m 2 
1 



\/ {2 ' 0) fa;) 

2 {p^,^^x^(x.<?E)} 



16m 

i vSZHx 



X, 



^{Pr,xx fl B}, 

x 



t^ o^ = ot^r o^r + otc 2 ' o) o«' o) (Px), 

O^Oh^ = —V^x) OV 1 ) . 0gB 



with 



1 ^OgiOcfo) 



._^M t(x.„")0(x. 9 B), 
O'COC'IP.) = !^<S!«M [0V> ■ (P_vOx, 9 E)] 

(x.^JPAfjExx)} 



(33) 
(34) 

(35) 



16m 2 


X 2 


1 




16m 2 


X 2 


1 




16m 2 


X 2 



{O f P x • (x x crW) 



crW)0,(x.^E)} 
^{Ot(x.P>WO,(xx 9 E)}. 



(36) 

We have defined Px = -«Dj and p x = — iV x = (pi — P2)/2; when acting on a 
singlet field Px reduces to -iVx = pi + P2 5 Pj being — iV x .. The abbreviation C.C. 
stands for charge conjugation, H.C. for Hermitian conjugation. For the potential 
we have used a notation close to the traditional one of [14] : we indicate momentum- and 
spin-independent potentials with V r , momentum-dependent potentials with V p 2 and Vl,2, 
spin-orbit potentials with Vls, spin-spin potentials with V52 and spin-tensor potentials 
with Vs 12 . The letters a, b, ... that appear in some of the matching coefficients are used 
to label the different kinds of operators. For ft/ 1 ' 1 ) and h^ 2 ' ^ only the Px-dependent 
terms are displayed. 

Under color trace all displayed terms of the type hg'p vanish. Due to charge conjuga- 
tion the terms proportional to V^, V<$$>, V^J , V^fJ , V^fJ,, V^fJ,, and V^fJ vanish. 
The matching coefficients c^~ 2 \ c^°\ V,i°/\ V%$, V&°] , and V { ^ ] are 



OOa i v OOa i v tp<f>b i v (p<p 



<p4>a 
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equal to 1 at tree level. All other matching coefficients are zero at 0(a®). The tree-level 
matching has been performed by multipole expanding the NRQCD Lagrangian (10) and 
projecting on singlet and octet two-particles states. 



3.1 Canonical quantization 

The canonical variables and their conjugates are (S,iS^), (O a ,iOl), and (A ia , II* = 
d£ P NRQCD/<9 (^oA")). The physical states |phys) are constrained to satisfy the Gauss 
law: 

(D.nnphys) = J d 3 xTr{0^gT a , O]} |phys). (37) 
The canonical commutation relations are 

[Ui(X,t),Ai(Y,t)\ = iS^jS^iX-Y), (38) 
[S(x,X,t),S\y,Y,t)} = «f< 3 >(x-y)<f< 3 >(X-Y), (39) 
[O a (x,X,t),Ot(y,Y,t)} = 5 ab 5 (3 \ X -y)5^(X-Y), (40) 

all the other commutators are zero. As in the case of NRQCD, in order to fulfil the 
Gauss-law constraint (37), we will assume the Poincare algebra commutators to act on 
a space spanned by the physical states. 



3.2 Poincare algebra generators in pNRQCD 

As in the case of NRQCD, since translational and rotational invariance are exact sym- 
metries of the effective theory, the generators H, P and J can be derived from the 
symmetric energy-momentum tensor. The pNRQCD Lorentz-boost generators K can be 
derived by writing down the most general expression, with some obvious restrictions al- 
ready included, containing all operators consistent with its symmetries and by matching 
it to the NRQCD Lorentz-boost generator (20). 

In our case we obtain at order x 2 /m°, x°/m, (x/m) Px and (x /m 2 ) P x (and drop- 
ping terms involving 4 matter fields that show up at O(x 2 /m ), but shall not affect our 
results): 

11 B + f d 3 xTi (s t (2m + h s )S + O f (2m + h o )0 



h 



] h* Oh B 00 



+ [{$h so O + H.C.) + C.C.J + \p%ooO + C.C 
H = J d 3 X h. 

The generators P and J are exactly known: 

P = J d 3 X J d 3 a;Tr js+PxS + O+Pxoj + d 3 X [n a x,B a ], 

J = J d 3 X J d 3 xTi |s f x P x + x x p x + ^°" (2 ^ S 



C.C. 



E 



= n 

(41) 
(42) 



+O f X x P x + x x + 1 O 



d 3 XXx[n a x,B a ]. (43) 
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The Lorentz-boost generators at order x 2 /m°, x°/m and (x/m) Px are given by: 

K = -t P + J d 3 X ^ {X, /i} + y d 3 X J d 3 x Tr { [&k ss S + C.C." 

+ [(S^soO + H.C.) + C.C.] + [O+kooO + C.C.] }, 

W^-V+^+tih^Pr), with 
1 



1,(0,2)* _ 
1,(1.-1) i 



1 f H^) 



r 



8m 
1 

8m 
1 k {1 '~ 1] 

1 | K i4>4>b"' 

8m x 2 



^HK'xp,)'} 



1 ,- ^ ( x ) 

o,™ I 



r 



(x) 



+^{C ) "(-)-M 5 xpL} 



1 r ^ii^ ( X ) i 



'''Sm^ 



x £ x fe ,P^ p^;} 



ar 



8m 

8m x 2 

8m x 2 e ^ XXCT P ^ 

8m x 2 3 X) 



x 



C 1)l (Px) 



0, 



(44) 

(45) 
(46) 



(47) 



(48) 
(49) 



where kffl are matching coefficients specific of K. For k^ 1 ' 1 ) only the Px-dependent 
terms are displayed. 



The terms of K proportional to lc S g , k 



(0,2) ,,(1,-1) 
SO i 



and kg 1 ^ vanish under color trace. 



Due to charge conjugation, also the term proportional to k^ 2 •* vanishes. The matching 
coefficients &ooa> fc™ 1 ') ^i^a' an d are equal to 1 at tree level. All other matching 
coefficients, which are specific of K, are zero at 0(a®). Also in the case of K the tree- 
level matching has been performed by multipole expanding the NRQCD Lorentz-boost 
generators (20) and projecting on singlet and octet two-particles states. We note that 
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loop corrections can , in principle, be calculated as they are for the matching coefficients 
of the pNRQCD Lagrangian. 



3.3 Poincare algebra constraints in pNRQCD 

Let us now consider the constraints induced by the Poincare algebra (l)-(9) on the 
pNRQCD generators H and K. 

As in the NRQCD case, the constraint [P l , K- 7 ] = —i5ijH has been already used 

in writing Eq. (44). Indeed, it forces K to have the form J d 3 X {X, h(X, t)}/2 + 

translational-invariant terms that depend on X only through the canonical variables. 
From [K\ K*] = -ie ijk J k at O(x /m°) it follows that 

, (1,0) _ , (1,0) _ , (1,0) _ , (1,0) _ 1 ,, n x 

K SSa' ^SSa" ~ K OOa' K OOa" ~ x > l ou I 

, (1,0) _ , (1,0) _ 1 r , n 

K ssc — K ooc — L i i 01 ) 

, (1,0) _ , (1,0) _ n r , 2 x 

^SSd' — K OOd' — u ; \ oz ) 

, (1,0) , (1,0) _ , (1,0) , (1,0) _ n ( ro\ 

K SSd" "T K SSd"' — K OOd" "T K OOd"' — u > \ 00 ) 

cg' 0) = cg< 0) = 1. (54) 

The constraint (54) follows also from [H, K. 1 ] = —iP % at O(x /m°). 

From [K\ K*] = -ie ijk J k and [H, K?] = -iP* at O(x/m ) it follows that 

= (55) 

= Vfe". (56) 

From = -iP* at O(x /m), we obtain: 

y (l,0) _ y (2,0) /roN 



VL5 5a + (fcg^ + = 0, (59) 

1/(0)/ 

VLsro. + + = 0, (60) 

t/(0)' 

Vfc o a + (*» + kg& + kSSl) X \ = 0, (62) 



T/(0) T/(0)/ 

Vp 2 Sa + VL**. + ^- - *ffi„ ^f - = 0, (63) 
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t/(°) , , r v<®' 

Vp 2 Oa + Vvoa + -f- ~ fcg&'-f- = 0, (64) 

ACFVsOb °s v SOa K SSa V SO ~ U > l D0 J 

7 (1,-1) / (1,0) 

^CfV/ 005 - c s v 00a V 00 — - U, (00 J 

7 (1,-1) 7 (1,0) 

9 T/( 1,0 ) T/( 2 '°) ft OOa T/(0,1) , ft OOc _ n f fi7 A 

A v o®ob - v otg>oa 2 — 00 ~^ 2~ ' 

- 4V." = 0, (68) 



2r£?-^-*&W = 0, (69) 

, (1,-1) T/ (0,1) , (1,0) 

9 T/^' ) T/( 2 '°) ^OOfe' K oo ft OO^ _ n / 7n \ 

z ^PPc ~~ ^PPfc' o n — — U > V U J 



^(1,-1)t/(°,1) 7.(1,0) 

>b' 

,(1,-1) 7,(1,-1) 



9 T/( 1 ' ) T/( 2 -°) ^OOb' ^OO , K OOd> _ n / 71 x 



k K ~ > - k ( ~ ' — (79) 

^SSb' K OOb' — u ; I ' Z J 

t/2,0) , , (1,-1) t/0,1)/ , (1,-1) / T/(0,1)Y _ n / 7 o\ 

^SOfe" + %5a xv SO + K SSb" \ XV SO ) ~ U > I'^J 

,(1,-1) T/ (0,l)/ 1.(1 -1) ( T T/^'^V 7(1,0) 

^oofc" 2 2 2 — ~~ 

,(i,-i) T /(o,i)/ t' 1 " 1 ' fTV (0,1 'V , (i,o) 

^p^pfe" 2 2 2 ' 

kssb"^ ~ ^obb" = 0) (76) 

^(2,0) , (1,-1)^(0,1) _ r?7 x 

v SOb"' + ^SSft'" V 5P — u ; I ' ' J 

, (1.-1W0.1) ,(1,0) 

T/ (2,0) . fc PPb"' ^PP K OOd"> n / 7S \ 



,(1,-1) T/ (0,1) u (lfi) 
y OOb 

,(1,-1) 7.(1-1) 



T/ (2,0) . ft PPb"' ^PP K OOd'" _ n / 7Q \ 

^o^ofe'" H 2^ 2 — ~~ 



^ssf/" — kobb'" — 0) (80) 

y(l,0) _ (1,-2), (0,2) , (1,0) „(2,0) _ n f81 x 

^PPa C P ft PPa + z ^PPa" + V PPc' — U ; ( 0i J 

^(1,0) (1,-2) , (0,2) _ , (1,0) _ v (2,0) _ n f82 x 
^PPa + C P K PPa Z K PPa' ^PPc" ~ U 5 

(1,-2) K OOb) (i,Q) _ y (2,0) _ n ,ooN 

Cq Z H 00 a'" v OOc"> — U ; (°°J 

r (l," 2 ) T i,(0,2)' _ 2 , (1,0) _ y(2,0) _ p , 84 x 

C G X K 00a ^^oOb V OOd — U 5 l S4 J 

where /' = df/dx. More precisely, as in the case of NRQCD discussed in Sec. 2.3, 
what we obtain are the above combinations of matching coefficients multiplying some 
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pNRQCD operators. For simplicity, here we have not displayed the expressions involving 
the operators. However, they may turn out to be important in considering the QED case 
of the above equalities. So, for instance, the expressions appearing in Eqs. (66) and (67) 
contribute to the same operator in the QED case. Therefore, the two equalities reduce 
to their sum, which is given by Eq. (65). Same considerations apply to the couples of 
equations (70)-(71), (74)-(75) and (78)-(79). Moreover, the combinations of matching 
coefficients appearing in Eqs. (81)-(84) multiply operators that vanish in the QED case. 
Therefore, these constraints have no QED analogue. For what concerns the equalities 
(55)-(64) in the QED case, they become pairwise identical. 



3.4 Unitary transformations 

The generators of the Poincare algebra are defined up to unitary transformations. We can 
use this freedom in order to change or reduce our basis of operators. (For a similar use 
of unitary transformations in a slightly different context we refer to [19].) In particular, 
we can look for a basis where some of the matching coefficients of K are fixed. Consider 
the unitary transformation: 



U 



exp 



(fX / d 3 xTr 



st 



\ K ss + K ss 



+ 



X 



2m 



s + c.c. 



2m 



o + c.c. 



(85) 



where the operators k^" 1 are equal to the operators but with arbitrary coefficients 

The transformed Lorentz-boost generators WKW^ have the same structure as the 
original ones up to order x 2 /m°, x°/m and (x/m) Px, but with the matching coefficients 
shifted as follows: 



1.(1,-1) 

1.(1-1) 
'></,&" 

1.(1.-1) 

1.(1,0) _ 
"V<K 

1.(1,0) . 
'V</>a" 

£,(1,0) 
1.(1,0) _ 

K tp<pb 

1.(1,0) 

1.(1,0) 
K <p<t>d"' 



> frt 1 " 1 ) 4- ~k (1 '- 1] 


(86) 


. 1.(1.-1) 7 (1,-1) 


(87) 


7(1,-1) r(i,-i) 
* K <fi4>b" K tp4>b" > 


(88) 


7(1,-1) r(i,-i) 


(89) 


, (1,0) , r (i,o) r (i,o) 


(90) 


, (1,0) r (1,0) r (1,0) 


(91) 


. 1.(1,0) , 2 r.(i,o) 


(92) 


i.(i.o) , 2 £(1.0) 


(93) 


, (1,0) r (1,0) r (i,o) 

"V0d" "Vi/id" "V^d'"' 


(94) 


, (1,0) 7,(1,0) 7,(1,0) 
^V^d'" "V0d" ^ip4>d'"- 


(95) 
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We have, now, the freedom to choose the coefficients . A convenient choice is that 
one that fixes the above matching coefficients to their tree level values: 

, (1,-1) , (1,0) 1 , (1,-1) , (1,-1) , (1,-1) , (1,0) , (1,0) , (1,0) , (i,o) , (i,o) n (q(]] 

The unitary transformation leaves invariant the generators P and J, while the changes 
in U HU^ may be reabsorbed in a redefinition of the matching coefficients. 

Finally, we note that other unitary transformations are possible and some of them 
may, in principle, further reduce the number of operators in K or H . We will not examine 
this issue here, which, however, deserves further studies. 



3.5 Discussion 

We comment now on the obtained relations. Equation (54) fixes the centre-of-mass 
kinetic energy to be equal to P^/4m. It may come as a surprise that, at least at the 
order at which we are working here, Poincare invariance does not fix the coefficient 
of the kinetic energy p 2 /m of the quarks in the centre-of-mass frame. However, one 
should consider that the coordinate x is no more a dynamical variable of the theory. 
Nevertheless, we may argue that, because no other momentum-dependent operator than 
the kinetic energy of NRQCD, — ^ V 2 / (2m) ip + x^ V 2 / (2m) x, may contribute to the 
kinetic energy of pNRQCD, the coefficients eg and Cg' 2 ^ have to be equal. It follows 
then that also Cg 1 ' -2 ^ = 1 (analogously for Cq'~ 2 ^). 
From Eqs. (59) and (60) we obtain: 

VLSSa _ VlSOa _ , / q7 x 
T/(0)/ _ T/ (0)/ - X ' ^'J 

where the last equality holds in the basis of operators discussed in Sec. 3.4. Eq. (97) in 
the singlet sector is the relation between the spin-orbit potentials and the static potential 
first derived in [10] by boosting the potentials expressed in terms of Wilson loops. It 
was also obtained in [9]. In all these previous cases the Lorentz-boost generators were 
written in the basis of Eq. (96). The extension to the octet sector is new. 
From Eqs. (61)-(64) we obtain: 

rV (o)> rI/ (o)' 
Vvsa + ~\- = Vvoa + ~\- = 0, (98) 

y(0) y(0) 

V p 2 Sa + V L 2 Sa + ^- = V p 2 Qa + V L 2 Sa + = 0. (99) 

These equations hold in the basis of Sec. 3.4. In the singlet sector they are the relations 
between the momentum-dependent potentials first derived in [11] by boosting the poten- 
tials expressed in terms of Wilson loops. They were also obtained in [9]. The extension 
to the octet sector is new. 

Equations (55)-(58) constrain the fields to enter in the singlet-octet and octet-octet 
sectors of the Lagrangian just in the combination 

x '( 9E+ Mlf x ' 9B })- < 100 » 
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i.e. like in the Lorentz force. The coefficients VgQ 1 ^ and Vqq^ were called Va and 
Vb respectively in the previous literature [20, 6]. These are the coefficients associated 
to the next-to-leading order terms of the multipole expansion of the static pNRQCD 
Lagrangian. They play an important role in the running of the static potentials [20, 21] 
and in several observables [18]. Eqs. (65)-(80) mix them with the matching coefficients of 
operators appearing at order 1/m and 1/m 2 . Eqs. (73)-(75) involve also the derivatives 
of VgQ 1 ^ and Vqq\ Eqs. (65) and (66) contain combinations of matching coefficients 
inherited from NRQCD. Somehow these relations reflect at the level of the pNRQCD 
potentials the relation (23) among the NRQCD matching coefficients. 

The four sets of Eqs. (65)-(68), (69)-(72), (73)-(76) and (77)-(80) can be combined 
in order to give the following four equations that involve only potentials appearing in 
the Lagrangian: 



T ,(l,0) T ,(2,0) o (rr,V {1 ' 0) 4 V {1 ' 0) \ - (r \/ (2 ' 0) 4- \/ (2 ' 0) ^ 
2c F V^Q b ' - C s V£o a > _ 1 \ c FV oob + V ®Ob) [ C s V OOa + V 0®Oa) 

T/(0,1) ~~ T/(0,1) 

v so v oo 
T/ (i,o) T/ (2,o) o ( V {1 ' 0) 4 1/ (1 ' 0) - ( V {2 ' 0) 4 V {m \ 

2 KsOc _ V SOb> _ \ °° c + V 0®Oc J { v OOb> + v O®Ob' J 

v so v oo 

V {2 ' 0) + ?ml/ (1 ' 0) - r I/ (2 ' 0) 
V SOb" + ZC F VsOb C s v SOa 



(101) 



0, (102) 



v (2,0) _ v (2,0) 9 ( yO-,0) , T/(1,0) \ / ^(2,0) ^(2,0) \ 

v OOb" v O®Ob"^ z \ c F v OOb + v O®Ob) \ c s v OOa + v O®Oa) 



(103) 



,(2,0) T/ (2,0) T/ (2,0) 



v SOb>" _ v OOb"> + v O®Ob'" _ n (-\(\A\ 

t/(o,i) - T/ ( ,i) ~ u ' 

v so v oo 

where the last equalities hold in the basis of Sec. (3.4). 

Eqs. (81)-(84) are typical for the non-Abelian structure of QCD and have been 
commented above. Here, we add that, summing Eq. (81) and (82), we obtain: 



, ^ T/ (2 ' 0) T/ (2 ' 0) 

Vgg> = 1 + W W - (105) 



4 Conclusions 

In this work we have shown how the implementation of the Poincare algebra provides 
a feasible and useful tool to constrain the dynamics of non-relativistic effective field 
theories of QCD. The method was used in the past, at a quantum-mechanical level, to 
constrain the relativistic dynamics of some phenomenological potentials. This is its first 
use in an effective field theory context. 

In Sec. 2 we have applied the method to the NRQCD/HQET Lagrangian up to order 
1/m and (V 2 V*)/m 2 in the commutators. In a simple way we have obtained some 
of the relations, Eqs. (21)-(24), derived formerly from reparametrization invariance. 
This shows that reparametrization invariance is, indeed, one way in which the Poincare 
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invariance of QCD manifests itself in the HQET (for a similar observation in a different 
context see [22]). An obvious extension of this work would be to study the constraints 
induced by Poincare invariance at higher order in 1/m. This study would be of interest, 
because at some point the four-fermion operators of NRQCD will start to play a role and 
new relations, specific of NRQCD, will show up. To our knowledge, relativistic invariance 
of NRQCD, in the sector where it does not reduce to the HQET, has never been explored. 
We believe that the framework used here is suitable also for this exploration. 

In Sec. 3 we have calculated the constraints induced in pNRQCD by the Poincare in- 
variance of QCD up to order x/m° and x°/m. We have obtained a set of new constraints, 
listed in Eqs. (50)- (84). These constraints involve, in general, mixings of matching co- 
efficients appearing in the pNRQCD Lagrangian and in the Lorentz-boost generators. 
There are two kinds of information that we can extract from them. First we can com- 
bine them in order to obtain relations that involve only the potentials appearing in the 
pNRQCD Lagrangian. These equations are (54)-(58), (97)-(99) and (101)-(105). More 
relations of this kind are expected in going to higher orders. In some specific bases 
of operators the relations become particularly simple. In the present letter we have 
presented a unitary transformation that fixes the matching coefficients specific of K at 
their tree-level value. An exhaustive analysis was beyond the purposes of the present 
letter but surely deserves further investigations. A second information, provided by the 
constraints, is the form of the Lorentz-boost generators expressed in terms of the poten- 
tials appearing in the Lagrangian of pNRQCD. Once K is known, the singlet and octet 
quark- ant iquark fields transform under infinitesimal Lorentz boosts with velocity v in 
accordance to 6S = i[S,\- K] and 50 a = i[O a , v ■ K]. 

Finally, we would like to mention that the present approach is quite general. There- 
fore, it may be also suited to derive exact relations among the matching coefficients of 
other effective field theories of QCD, where the manifest covariance under boosts has 
been destroyed by an expansion in some small momenta. An example is the soft-collinear 
effective theory [23], where also reparametrization invariance has been discussed [24]. 
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